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Abstract 

Gauge Higgs Unification in Warped Extra Dimensions provides an attractive solution to the 
hierarchy problem. The extension of the Standard Model gauge symmetry to SO (5) x U(l)x 
allows the incorporation of the custodial symmetry SU(2)r plus a Higgs boson doublet with 
the right quantum numbers under the gauge group. In the minimal model, the Higgs mass 
is in the range 110-150 GeV, while a light Kaluza Klein (KK) excitation of the top quark 
appears in the spectrum, providing agreement with precision electroweak measurements and a 
possible test of the model at a high luminosity LHC. The extension of the model to the lepton 
sector has several interesting features. We discuss the conditions necessary to obtain realistic 
charged lepton and neutrino masses. After the addition of an exchange symmetry in the bulk, 
we show that the odd neutrino KK modes provide a realistic dark matter candidate, with a 
mass of the order of 1 TeV, which will be probed by direct dark matter detection experiments 
in the near future. 
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1 Introduction 



Warped Extra Dimensions present an elegant solution to the hierarchy problem, where all funda- 
mental parameters are of the order of the Planck scale. The weak scale-Planck scale hierarchy is 
obtained by an exponential warp factor, which is naturally small provided the Higgs field is local- 
ized towards the so-called infrared brane pp. If all Standard Model (SM) fields propagate in the 
bulk, the theory leads to the presence of Kaluza Klein (KK) modes which tend to be localized 
towards the IR brane and therefore couple sizably to the Higgs. This in turn leads to large mixing 
between the heavy SM particles and their KK modes, leading to modifications of the electroweak 
parameters and therefore to strong constraints from electroweak precision measurements [2]-[7j. 
These constraints may be weakened by the introduction of brane kinetic terms [8] [11] or custodial 
symmetries [12], [13], which allow the presence of KK modes with masses of the order of a few TeV. 

One of the attractive features of these models is the natural explanation of the hierarchy of 
fermion masses by the localization of the fermion fields in the bulk [T6]-[T8]. The chiral properties 
of the fermions are obtained by imposing an orbifold symmetry and demanding that the fields are 
odd or even under such a symmetry. Fermion fields that are even under the orbifold symmetry 
at the infrared and ultraviolet branes present zero modes, which are chiral and therefore may be 
identified with the SM fermion fields. The localization of the zero modes is governed by the bulk 
mass parameter ck, with k the curvature of the extra dimension and c, a number of order one. While 
the zero modes of chiral fields with c > 1/2 couple weakly to the Higgs and are therefore light, the 
heavy SM fields are associated with bulk mass parameters c < 1/2. Due to the exponential behavior 
of the zero mode wave functions, large hierarchies between the fermion masses are generated by 
small variations of the corresponding c-parameters. 

Gauge Higgs Unification models identify the Higgs field with the five dimensional component of 
the gauge fields |I9] . An extended gauge symmetry is necessary for the successful implementation of 
this mechanism. In particular, models based on the gauge group SO(5) x U{l)x include the custodial 
and weak gauge symmetry via SO (5) D 50(4) = SU(2) L x SU(2) R [20] , [21] [22]. Moreover, 
provided the SO (5) symmetry is broken to a subgroup of 5*0(4) by boundary conditions at both 
branes, the fifth dimensional components of SO(5)/ SO(4) have the proper quantum numbers to be 
identified with the Higgs field, which is exponentially localized towards the IR brane. 

Since the Higgs originates from gauge fields, its tree level potential vanishes. In a previous 
work [25], we computed the one-loop effective potential and demonstrated that electroweak sym- 
metry breaking, with the proper generation of third generation quark and gauge boson masses may 
be obtained for the same values of the bulk mass parameters that lead to agreement with precision 
electroweak data at the one-loop level. Moreover, we showed that the Higgs mass is in the range 
110-150 GeV and that a light KK mode of the top quark, T", appears in the spectrum, with a mass 
small enough so that the KK gluon may decay into it. The presence of this light KK mode has a 
strong impact on the phenomenology of the model [26]. For instance, searches for the KK gluon 
by its decay into top-quarks [27], [28], [29J is rendered difficult due to the presence of the additional 
dominant decay mode into KK top-quarks and the broadness of the KK gluon. On the contrary, the 
constructive interference between the QCD and KK gluon induced pair production of the top-quark 
KK mode allows to search for a T' for values of the masses much larger than those at reach in the 
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case of just the QCD production cross section. 

In this article, we will analyze the addition of the lepton sector in the gauge Higgs unification 
scenario. To add leptons, we will proceed in a similar way as for the quark sector. The left-handed 
leptons will be added in a fundamental representation of 5*0(5), with Qx = 0, while the right- 
handed neutrino and charged lepton fields will be added in a fundamental representation and a 10 
of 50(5) also with charge Qx = 0, respectively. 

Due to the gauge origin of the Higgs field, a possible local infrared brane operator (LH) LH/M, 
which could lead to large values of the neutrino Majorana masses, should come from the fifth 
dimensional component of the covariant derivative of the lepton fields and therefore can only nat- 
urally arise from the integration of the right-handed neutrinos, with a local IR Majorana mass. 
Indeed, since the fields associated with the right-handed neutrino zero modes are singlets under the 
conserved gauge groups on both the infrared and ultraviolet branes, one can always add Majorana 
masses for these fields on the infrared and ultraviolet branes. We will therefore consider these masses 
and implement a See-Saw mechanism for the generation of the light neutrino masses [2], [32]. We 
will show how to incorporate these masses within the context of these models, obtain the modified 
profile functions and define the conditions to derive a realistic spectrum. 

Furthermore, the introduction of an exchange symmetry [13J, which is preserved in the bulk, 
yields a natural dark matter candidate in the spectrum, that may be identified with the odd KK 
modes of the right-handed neutrino. This exchange symmetry requires the identification of the 
bulk mass parameter of the dark matter candidate with the one of the right-handed neutrinos, 
establishing a connection between neutrino masses and the relic dark matter density. We will show 
that, if the odd fermions are assumed to be Dirac particles, the predicted relic density is the correct 
one for Dirac fermion masses of the order of 1 TeV. In the Majorana case, somewhat smaller masses 
are allowed, and the results depend on the relative values of the ultraviolet and infrared Majorana 
masses. We will show that in both the Dirac and Majorana cases, direct detection experiments will 
efficiently probe the existence of the proposed dark matter candidates. 

The article is organized as follow: In section 2 we discuss the properties of leptons in scenarios 
of Gauge Higgs Unification. In section 3 we discuss the generation of charged lepton and neutrino 
masses. In section 4 we analyze the possibility of incorporating dark matter via an exchange bulk 
symmetry, analyzing the couplings and the annihilation diagrams, as well as the direct detection of 
these dark matter candidates. We reserve section 5 for our conclusions. Some technical details are 
given in the Appendix. 

2 Leptons in Gauge Higgs Unification Scenarios 

The goal of the current work is to add the lepton sector into the minimal Gauge Higgs Unification 
model described above, including right handed neutrinos, with brane Majorana masses. We will 
proceed in a similar way as in the quark case [25], and discuss the general question of charged lepton 
and neutrino mass generation. 

Similar to the quark sector, we let the SM SU{2) L lepton doublet containing the left-handed 
charged leptons and neutrinos, II and z/ L , arise from a 5 of 50(5) x U(l)x, where the subscript 
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refers to the U(l)x charge. The right handed neutrino, N R will be included as the singlet component 
in a fundamental representation of SO (5), while the right-handed charged lepton, l R , is placed in a 
10o of 5*0(5) analogously to the d R . One may also include brane mass terms connecting different 
multiplet components, as well as new brane Majorana masses for the right handed neutrino N R : 

(M uv S(y) - M m 5(y - L))n r N r . (1) 

The right handed neutrino, N R , in principle could be identified with the singlet right-handed 
component in the same multiplet as the left-handed leptons. However, in order to naturally suppress 
lepton flavor violation effects and maintain agreement with precision electroweak measurements, 
the left-handed leptons should be localized towards the UV brane [I4].|15|. The zero modes of the 
corresponding right-handed multiplet will therefore be localized towards the IR brane. This implies 
that even with a natural scale for the brane masses O(Mpi), the exponential suppression of the wave 
function at the UV brane would lead to a effective Majorana mass for N R which is much smaller 
than the Planck scale. This, in turn, after the implementation of the See-Saw mechanism, leads 
to too large values for the observed neutrino masses. Therefore, as we will discuss below, it will 
prove to be necessary to have the left-handed leptons in a different multiplet as the right-handed 
neutrinos. 

If N R belongs to the same multiplet as the left-handed leptons: 

^ LlL ~ Ui(-,+)o l e £(+,+U ® iY ^-'-)°' (2) 
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Alternatively, the right-handed neutrino can be incorporated in a different multiplet from the 
left-handed lepton one. The two multiplet assignments £1X6 cLS follows: 
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where we show the decomposition under SU (2) l x SU(2)r, and explicitly write the U{1)em charges. 
The Us are bidoublets of SU{2) L x SU{2) R , with SU{2) L acting vertically and SU{2) R acting 
horizontally. The T"s and T 2 s transform as (3, 1) and (1, 3) under SU{2)l x SU(2)r, respectively, 
while N % and N' % are SU(2)l x SU{2)r singlets. The superscripts, i = 1,2,3, label the three 
generations. 

We also show the parities on the indicated 4D chirality, where — and + stands for odd and even 
parity conditions and the first and second entries in the bracket correspond to the parities in the UV 
and IR branes respectively. Let us stress that while odd parity is equivalent to a Dirichlet boundary 
condition, the even parity is a linear combination of Neumann and Dirichlet boundary conditions, 
that is determined via the fermion bulk equations of motion as discussed below. The boundary 
conditions for the opposite chirality fermion multiplet can be read off the ones above by a flip in 
both chirality and boundary condition, for example (— , +)l — > (+, —)r. In the absence of mixing 
among multiplets satisfying different boundary conditions, the SM fermions arise as the zero-modes 
of the fields obeying (+, +) boundary conditions. The remaining boundary conditions are chosen 
so that SU(2) l x SU {2)r is preserved on the IR brane and so that mass mixing terms, necessary 
to obtain the SM fermion masses after EW symmetry breaking, can be written on the IR brane. 
Consistency of the above parity assignments with the original orbifold Z 2 symmetry at the IR brane 
was discussed in Appendix B of Ref . [25] . The three families will behave similarly, and therefore, we 
will drop the family indices and concentrate only on one lepton family. Large mixing angles in the 
lepton sector can be naturally obtained while suppressing lepton flavor changing neutral currents 
if the left-handed leptons have similar bulk mass parameters, c\ [36], [37], and in the following we 
will assume them to be equal. We will return to this issue in section 13.31 The zero modes of the 
leptons are too light and too weakly coupled to the Higgs boson to affect the Higgs potential in any 
significant way. The lepton KK modes may be coupled more strongly to the Higgs, but their gauge 
invariant mass is much larger than the Higgs induced one and therefore they tend to contribute 
only weakly to both the Higgs potential and to precision electroweak observables. 

One can add localized brane mass terms to the Lagrangian in both the one and the two multiplet 
cases: 



C 2 



-25(x 5 - L) 



L 1L M L2 L 2R + h.c. 



- [M IR 5(x 5 - L) - M uv 6(x 5 ) 



N R N R ; 



(4) 



25(x 5 - L) N' L M Ll N R + L 1L M L2 L 3R + h.c. - [M IR 5(x 5 - L) - M uv 6(x 5 ) N R N R (5) 



With the introduction of the brane mixing terms, the different multiplets are now related via 
the equations of motion. The fermions, like the gauge bosons, can be expanded in their KK basis: 



1p L (x, X 5 ) = ^ fL,n( X 5, h)^L,n{x), ^r{x, X 5 ) = ^ fnA X 5, h)lp R) r 



X 



(6) 



Solving the equations of motion in the presence of h becomes complicated, as the different modes 
are mixed. However, 5D gauge symmetry relates these solutions to solutions with h = [23], with 
f2(#5, h), the gauge transformations that removes the vev of h: 

rx 5 



fl(x 5 , h) = exp 



-iC h hT 4 



dya (y) 



(7) 
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The form of the bulk profile functions at h = is given in Appendix A. 

The boundary masses lead to a redefinition of the effective boundary conditions for the fermion 
fields at the branes. Let's first analyze the case of a Dirac boundary mass term on the infrared brane, 
Ml, involving fields from different multiplets ty l L ty J R , i, j = 1, 2, 3 for the different multiplets. Quite 
generally, we shall call the profile of the left-handed field and the right-handed field participating in 
the mass term and hji, respectively. The right handed component, ty R with profile function g R 
and the left-handed component *ff 3 L with profile function hi have Dirichlet boundary conditions on 
the brane, and therefore g R (L) = hi(L) = 0. The equations of motion of the fields are affected by 
the localized masses, which induce a discontinuity on the odd-parity profile functions at the infrared 
brane. Indeed, keeping only the relevant terms, the integration of the equation of motion leads to 

/ (d 5 g R )dx 5 = / 2M L h R 6(x 5 - L)dx 5 . (8) 

JL-e JL-e 



Therefore, we obtain: 



and, similarly for the 



\img R (L-e) = -M L h R (L), (9) 

e— >0 



\imh L {L-e)=M L2 g L {L). (10) 

Eq. (J9j) and ( jTUl) can now be reinterpreted as the new boundary conditions for the profiles at the 
IR brane. 

Analogously, one can analyze the effect of the Majorana boundary mass, M», where % = IR or 
UV. Let's take the specific case of the field N R , with a profile function h R . Its chiral partner will 
have a profile function hi having an odd parity profile on both branes. The equation of motion 
in the presence of both the Majorana masses and the Dirac mass term Ml 1 leads to the following 
relationship 

[ V+ (d 5 h L )dx 5 = [ V+ U±2Mih R - h L )5(x 5 -y)- 2M Ll g L 5(x 5 - L)j dx 5 . (11) 

where the minus and plus signs are associated with the boundary conditions at y = L, and y = 
0, respectively. The odd parity at the branes then implies a Dirichlet boundary condition for 
the function hi, which as before will present a discontinuity at the brane. For the IR boundary 
conditions we obtain: 

\imh L (L - e) = M IR h R (L) + M Ll g L (L). (12) 
For the UV boundary condition, instead: 

hmh L (Q + e) = M uv h R (Q). (13) 

Eq. (I12H13I) can now be reinterpreted as the new boundary condition for the profiles at the branes. 
The generalization of these expressions to the general case is straightforward. 
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2.1 Wave Functions in the Presence of UV Majorana Masses 

The wave functions defined in Appendix A, Sm and S-m are associated with Dirichlet boundary 
conditions at the ultraviolet brane for the left-handed and right-handed fields, respectively. In 
the presence of ultraviolet Majorana masses, however, the boundary conditions for the singlet 
component of the fundamental multiplets of 5*0(5) read 

/^(0, 0) = M uv fl R (0, 0) (14) 

where % = 1, 2 refers to the particular multiplet under analysis. In order to derive the profile function 
of these fields, let us first redefine the functions a 3 ^ 2 fL,R — > Jl,r where a = exp(— kx^). With this 
redefinition, these functions satisfy the naive normalization condition, 

L 

d%5 fnfm = $rn,n- (15) 







The general solution for fi is given by: 

f L (x„ 0) = Aa^S M + B (~) a-^S'_ M - (16) 

where z is the associated particle mass. Defining 

h,R = a- (c " 1/2) /L,«, (17) 
fi t ji satisfy the simple equation of motion, 

Wi5,0) = T-Wi 6 ,l)). (18) 

Now, using Eq. (fl8|) . one obtains 

f R (x 5 ,0) = -^(A^-fi(^)a- 2c ((l-2c)^ M -^ M )). (19) 

The second derivative functions may be replaced by means of the equation of motion of the fermion 
fields, namely 

S1 M = k(l - 2c)S'_ M - Z -~S. M . (20) 



a 



We therefore see that fn reduces to: 



a ~ 



f R (x 5 , 0) = Ba~ 2c S^ M - A-S' M . (21) 



z 



Rewriting these in terms of the S rather than S', with S±m = T^^-S' ±M , we obtain 

f L (x 5 ,0) = AS M + Ba- 2c L M (22) 
f R (x 5 ,0) = Ba- 2c S„ M + AS M (23) 
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To solve for A in terms of B, we need to use the UV boundary condition induced by the UV 
Majorana mass for Nr, Eq. (fl4|) : 



Ba~ 2c S_ M (0) 



since a = l and S± M (0,z)=z for X5 



.4 

0: 
.4 



AM uv S M (0) 



B- 



-M 



M 



uv S M 
B 



1^5=0 



Mi 



(24) 



Therefore, with the coefficients A and B as calculated above, the singlet functions become 

fl L {x B ,0) = C 5 (S Ml -M uv S. Ml ) (25) 



fl R (x s ,0) = C 5 (-M uv S„ Ml + S Ml ) 
in the case of a single multiplet containing the left- and right-handed neutrinos, and 



/ 2 5 l(^5,0) 



C\q(Sm 2 — MuvS-m 2 ) 
Ci (—M uv S^m 2 + Sm 2/ 



(26) 



(27) 
(28) 



in the case of two multiplets, where S±m = cS- ±c 1 ' 2 'S±m and S±m = cS ±c 1 ^ 2 * > S±m- Therefore, the 
fermion multiplets with h = take the form 



/i,l(z5,0) 



CiSm! 

C2Sm 1 

C4S-M1 

Jl,L 



/3,fl(af5,0) 



C12S-MS 
C13S-M3 

CyjS-M 3 
CigS-M 3 
CxgS-M 3 



(29) 
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in the case of a single multiplet containing the left-handed and right-handed neutrinos, and 



/i,l(>5,0) 



/2,Z(Z5,0) 



CiSm! 
CzSmx 



C7S-M2 
C9S-M2 

/2,L 



fi 





-M 3 
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-M 3 
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-Ms 
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-Ms 




-Ms 


CiqS. 


-M 3 


C17S- 


-Ms 


CigS- 


-M 3 


C19S- 


-Ms 


C20S 


Ms 



(30) 



in the two multiplet case, where the Cj are normalization constants. 



3 Lepton Spectrum 

Applying the boundary conditions at £5 = L, taking into account the mass mixing terms from 
Eqs. (BJ and <$J5§ and using the procedure defined in Eqs. (l9~i)- ffT2l . we derive the conditions on the 
lepton wave functions f(L, h) in the presence of the Higgs field. In the case of only one multiplet 
containing both the left-handed and right-handed neutrinos one gets the following conditions at the 
IR brane: 



fi,L ~ M IR f* = 



4r 10 = 



f} : ,r ■ M,j\ur 4 

/fe" 4 - M L2 fl :r A = 
In the two multiplets, instead, one obtains: 

fi,n A + M l JI :r a = ./?.„ + M l JIr = 



(31) 



rl,...,4 
J3.L 



MlJIt 4 = fL - M l JIl - M IR ff 



2,R 







fx, 



2,L 

5,. ..,10 





= 



(32) 



where the superscripts denote the vector components. 

This defines a system of linear equations for the normalization constants Cj. Asking that the 
determinant of the functional coefficients of this system vanishes in order to get a non-trivial solu- 
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tion [23], one obtains the following relations: 



(33) 



2S 



Mi 



M^SmxS-Mz + SmxS-Ms — 



M^S-m^S-Mx + S^m 1 S-m 3 — Ml 2 S-M! sin 



Xh 

fh 



(34) 
(35) 



—M1 2 Sm 1 S 2 l M3 S-m^Sm! — e 2cikL MuvS-M 1 ) + Mz\r(1 — Sm^-M! + e 2cikL MuvS-M 1 S-M 1 , 



-M?J. 



S. 



M 3 



-S-. 



Sm^Sm! — MirSmi) + e 2cifci M[/v(l — SmiS-mi + ^irS-Mi^. 



Mi, 



s 2 



M 3 



-(Ml 2 M IR S„ M3 + S. 



M 3/ 



M 2 L2 S- M3 



_l_ e 2c ^ kL MuvS-M 1 + Sm^Sm-lS-Mx — SmxS-M! 
in the case of a singlet neutrino multiplet, and 

£ 3 



&Mi + e 2ci kL M uv S*_ Mi (Sm x S- Mi — SmxS-Mi, 

.2 

S-M 3 



sm 



Aft. 







-M 2 







S 5 -M 3 = 



M L2 Sm x S-M3 + Sm\S-Mz 



n 2 



2S 



M 3 



Mi^S-MsS-Mi + S-M1S-M3 



M 2 2 S_ Ml sin 



Xh 

fh 



(36) 

(37) 
(38) 
(39) 
(40) 
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M 3 



1 — Sm 1 S-m 1 )(Sm 2 ~ e2c2kL MuvS-M 2 )S^M 2 + M IR {1 — SmxS-Mx, 



(1 — Sm 2 S-m 2 + e2c2kL MuvS-Ai 2 S-M 2 ) + ^| 1 £mi'S'-Mi(1 ~ Sm 2 S-m 2 + e 2c2kL Muv S -m 2 S - 



M 2 , 



Aft 



— Sm 2 S-m 2 + e 2c2kL MuvS-M 2 S- 



M 2/ 



S 



Aft 



S-m 2 {Sm 2 — e 2 ° 2kL MuvS-M 2 ) 



+M IR (1 — Sm 2 S-m 2 + e2c2kL MuvS-M 2 S-M 2 , 



5_ 
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+ 



M 3 



Aft 



+Mtr(— 3 + 2Sm 1 S-m 1 + Sm 2 S~m 2 — e 2c2hL MuyS-M 2 S- 
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Aft, 



2M ir S^m 1 Sm 1 — + 2Smj_S-mi ~ Sm 2 S~m 2 + ^ 2c2kL M uv S^M 2 S-M 2 ] 
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Aft 



sm 
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M 2 M ir S_ Mr + M 2 S. 



sin 



Aft; 

ft; 



-i 4 







(41) 



in the case of two multiplets. In the above, for simplicity, we did not write the dependence on L 
and z and furthermore, we have used the Crowian: 



Sm{x 5 , z)S- M (x 5 , z) + S M (x5, z)S^ m (x 5 , z) 



1 . 



(42) 



The roots of the above equations define the values of z corresponding to the masses of the lepton 
zero modes and KK modes in the presence of the Higgs fields. 



3.1 Charged Lepton Masses 

Since the charged lepton masses are given by the mixing of the first and third multiplet via Ml 2 , 
the expression determining its mass is formally the same for the case in which the right-handed 
neutrino is in the same multiplet as the left-handed one as for the two neutrino multiplet case 
(Eqs. fl35|) and (|40p ). Additionally, since the lepton masses are much smaller then k, one can use an 
expansion of Sm for small values of z/k. As we shall show, the approximate functions we derive in 
this way agree very well with the full numerical investigation we carried out. We shall concentrate 
on values of c± £ 0.5, which are preferred to cancel flavor violation effects in a natural way and 
provide agreement with precision electroweak data [T5] . |36] . 

At small values of z compared to k, one can express the function Sm m the form 



i M ~ z I a-\y)e- 2My dy + 0(z 3 



(43) 



Using this in Eq. ( |35l) . we can solve for the mass: 



M L2 e^ kL Mf}J2 (I- Cl )( c §-I 



2( Cl -l)kL _ x) _ M 2^ Ci _ !)( e 2( C3 -i)fc _ Vj 



c 3 )(e 



(44) 



( e 2(±+c 3 )kL 
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If c\ > 0.5 and —0.5 < C3 < ci, this reduces to: 



U\ =M L2 e^ kL si 
For ci > 0.5 and c 3 > Cj, instead, 



f 1 - 



sin 



Finally, for the case ci > 0.5 and C3 < —0.5: 



A/i 

fh 



2(ci-^)(ca + i) 



2(c| 



1, 

4' 



(45) 



(46) 



M L2 e (1+C3 - Cl)fcL sin 



A/i 



[~ -ci)(c 3 + -) 



(47) 



where we have assumed that M L2 7^ 0. We note that in the above, the lepton masses depend at 
most linearly on M^ 2 . As shown in Fig. [TJ the above relations are verified by our numerical results. 
Realistic lepton masses may be obtained for e.g. for c\ ~ 0.6 and C3 ~ —0.55, —0.65 and —0.8 for the 
case of the tau, muon and electron, respectively. If a common value of c\ for the three generations 
is demanded, as explained above, and for values of Ml 2 of order one, as chosen in Fig. [fl the value 
of Ci is restricted to be in the range 0.5 < c\ 1z 0.75. Larger values of c\ become incompatible with 
the heavier charged lepton masses. 



3.2 Neutrino Masses 

The Neutrino masses are analyzed in a similar manner. First we look at the case in which the 
left-handed and right-handed neutrinos belong to the same muliplet case. For c\ > 0.5, 

( c _ i) e 2(i- cl )fci sin rMi 2 
' - : (48) 



M IR 

From Eq. (T4H1) we see that values of C\ ~ 1 would be necessary in order to get the correct values 
for the neutrino masses. However, values of C\ ~ 1 are strongly constrained in order to reproduce 
the proper charged lepton masses. Indeed, as we emphasized at the end of last section, the proper 
values of r and fi masses may not be obtained for c\ ~ 1. Therefore, we conclude that if all ci's 
are about the same, as preferred to obtain large flavor mixing naturally without inducing large 
lepton flavor changing effects [36] , two multiplets are required in order to obtain the correct lepton 
spectrum. 

In the two multiplet case, the dependence of the neutrino masses on the mixing with the third 
multiplet through M L2 is always exponentially suppressed. Therefore, we shall set M L2 = in the 
following approximate expressions. The approximate mass expressions, for 

• Ci > 0.5 and c 2 > l/(k L): 
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Figure 1: Region of ci, C3 parameter space consistent with the known charged lepton masses. The bands correspond 
to variations of the values of the k and Ml 2 parameters in the range 1.5 TeV ^ k < 5.5 TeV and 0.1 < Ml 2 ~ 5. 
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(49) 



ci > 0.5 and c 2 < -l/(k L): 



M lx( c i - \)e 2 ^- c ^ +c ^ kL sin 



Xh 
fh 



1 2 



Mi 



uv 



(50) 



C\ > 0.5 and c 2 ~ 0: 
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where in Eq. ([ST]) we have assumed Mjjy ^ —Mm. 

In the linear regime, (Xh/fh) 2 <C 1, these neutrino masses become proportional to the square 
of the Higgs vev, and show the characteristic See-Saw behavior governed by the brane Majorana 
masses. From the above expressions we see that it will only be possible to generate the correct order 
of the neutrino masses when c 2 > —0.4. Moreover, for c 2 > 0, the values of C\ are such that the 
correct heavy lepton masses cannot be generated. These conclusions are verified in our numerical 
work. We present the relevant parameter space in the c\ — c 2 plane leading to the correct order 
of the neutrino masses in Fig. [2j The width of the bands for the different masses corresponds to 
varying k and the different brane masses in the range indicated in Fig. [2J As indicated by the above 
expressions, we were not able to numerically find any solutions for c 2 < —0.5. Finally, although 
positive values of C2 are not represented in Fig. [2j the neutrino masses become independent of C2 
for C2 > 0, and therefore the values of c\ are the same as for C2 = 0. 

3.3 Flavor Problem 

In the above, we have not discussed the problem of flavor. It is well known that, if the effective 
Yukawa couplings have an anarchic structure, large flavor violating effects are induced, which may 
only be suppressed by pushing the KK masses to values above 10 TeV, excluding any possible 
phenomenology of warped extra dimensional models at the LHC, as well as any possible dark 
matter candidate coming from the KK modes (see, for example Ref. [32] as well as Ref. [31] for an 
alternative approach to this question). The problem stems in part from the fact that the Yukawa 
couplings and the bulk mass parameters are not diagonalized in the same basis, and therefore 
the quark mass eigenstates have flavor violating couplings with the gluon KK modes. A possible 
solution to this problem is to demand an alignment between the bulk mass parameters and the 
Yukawa couplings, as has been proposed in Ref. [35]. Flavor violation in the lepton sector can also 
be suppressed by a similar alignment [36J,[37J. This is equivalent to demanding that the bulk mass 
parameters obey the following relationships: 

c 3 = I + a^Yt 
c 2 = I + a 2 k 2 YjY u 

Cl = I + atfYtf + a„k 2 Y v Y2; (52) 

where Yi and Y u are the effective charged lepton and neutrino Yukawa couplings and ai, a u , a 2 
and a% are numerical constants and the q are now matrices where the off-diagonal terms mix the 
different generations. 

The Gauge-Higgs unification structure introduced above demands a redefinition of the above 
equations, since no explicit Yukawa coupling has been written. As can be seen from Eqs. (j4"5]) - (j4"T|) . 
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Figure 2: Region of ci, C2 parameter space consistent with the neutrino masses of interest: c\ > 0.5 and —0.5 < 
C2 < 0. The bands correspond to variations of the values of the parameters k, and Mjjvjr in the range 

1.5 TeV < k < 5 TeV, 0.1 < M Ll < 1.5 and 0.5 < M m , uv < 2.5. 

and Eqs. fj49l)-( !5T|) . the role of the Yukawa coupling is now being played by the boundary masses 
Ml 1 and Ml 2 . Hence, the above equations must be rewritten as 

c 3 = I + a 32 Ml 2 M L2 
c 2 = I + a 2X Ml i M Ll 

ci = I + a 12 M L2 Ml 2 + anM^Ml; (53) 

If <2i2 ^> an, the charged lepton masses would be diagonalized in the same basis as the bulk 
mass parameters, inducing minimal flavor violation in the lepton sector. In this case, all flavor 
violation will be associated with the charged currents, leading to values of the lepton flavor violation 
processes consistent with experiment for KK masses as low as a few TeV. As emphasized above, large 
mixing angles may naturally arise within this framework, if all left-handed zero modes localization 
parameters take equal values, namely when a 12 ~ an ~ |36j . 
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The results of Refs. [36], [37] show that it is possible to impose a flavor symmetry in this class 
of models such that no large flavor violation occurs. In this work, we will assume that this, or 
a similar [38] , [39] , [40] flavor protection exists, and will postpone a more detailed analysis of this 
question for future work. 

4 Dark Matter 

Dark Matter in warped extra dimensions was first introduced in Ref. [H] within a framework 
which solves the proton stability problem in unification scenarios. The introduction of a KK parity 
in warped extra dimensions, leading to a stable dark matter candidate, was further proposed in 
Ref. [12]. In this work, we shall proceed in a different way: Following Ref. [43J, we shall introduce 
an additional exchange Z 2 symmetry under which all the lepton multiplets introduced so far would 
be even. One can then define extra fermion multiplets, that will be chosen as the "odd" partners 
of the lepton multiplets. If this symmetry is preserved, the lightest odd particle (LOP) will be 
stable, and therefore can be considered as a possible dark matter candidate. In the framework of 
Ref. [43] the equality of the even and odd mass parameters was enforced by giving the original 
fermions, whose even and odd combinations form the even and odd fields, different charges under 
an extended U(l)x 1 x U(l)x 2 gauge symmetry. Since in our case the leptons are neutral under 
U(l)x this property does not hold. Additionally, contrary to Ref. [43J, we shall assume that the 
quark and gauge boson multiplets do not have odd partners. 

Even though, the structure of our model does not require the equality of the bulk masses for the 
odd and even fields, for simplicity, we shall assume that the bulk mass parameters are identified with 
each other and are controlled by the requirement of obtaining the correct small neutrino masses via 
the see-saw mechanism. Our assumption is equivalent to requiring that there are no off-diagonal 
bulk mass parameters mixing the original fields for which the Z2 exchange symmetry holds. 

In order to explore this possibility, we shall identify the multiplet containing the dark matter 
candidate with the odd partners of the second lepton multiplet containing the right-handed neutri- 
nos. As has been shown in the previous section, this demands values of C2 < 0, and therefore we 
shall require the bulk mass parameter of the dark matter candidate to be in this range. 

The exchange symmetry, introduced in Ref. [42] allows arbitrary boundary masses between even 
fields, necessary for obtaining the proper lepton masses, as well as between the odd fields. Boundary 
masses mixing odd and even fields are, instead, forbidden by this symmetry. On the other hand, 
the boundary conditions for the even and odd fields are independent of each other. Therefore, the 
main link between even and odd fields is through the identification of the bulk mass parameters. 
For simplicity, we shall choose the boundary conditions of the odd partners of the chiral first and 
third lepton multiplets, containing the left-handed and right-handed charged leptons, to be the 
same as the one of the even biodoublet components, (- +) and (+ -), respectively. For values of the 
bulk masses c\ > 0.5 and C3 < —0.5, these fields will be relatively heavy, with masses about a few 
times k, and decoupled from the Higgs. 

For the second multiplet, the odd fields, denoted by £ D will have the same boundary conditions 
as £2 for the even bidoublets. However, since the fifth component does not contain a zero mode, 
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it must have different boundary conditions on the IR and UV branes. That leaves us with two 
options for the left handed component of this field: (+, — ) or (— , +). The goal here is to consider 
the possibility of a neutral odd lepton, mainly singlet under the SU(2)i x SU(2)r symmetry, as a 
dark matter candidate. The singlet and the doublet states mix via their interactions with the Higgs 
field, which will act as a small perturbation to their masses. In order to make the coupling to the 
Higgs effective and to split the doublet and singlet masses, we shall choose the singlet right-handed 
field to have the same boundary conditions on the IR brane as the bidoublet left-handed field for 
at least one of the three odd partners of the second multiplets. Therefore, the boundary conditions 
for the odd multiplet containing the LOP are chosen to be: 

a ~ L « = c%(-',+)-i) 9 Ar « (+ -- ) °- (54) 

Regarding the other two generations of odd partners of the second multiplets, for simplicity, we 
will choose their singlet states to have the opposite boundary conditions from the one presented in 
Eq. (1541) . This would force their masses to be heavy for C2 < 0, ensuring that the multiplet with 
the boundary conditions given by Eq. (J541) . £ , would generate the LOP. In addition, small Dirac 
boundary masses may be included, which would allow a small mixing between the odd multiplets 
inducing decays of the heavier generation odd states to the LOP, through the weak gauge bosons 
and the Higgs boson. Even in the case of a very small mixing, due to the large mass differences, 
the lifetime of these heavy odd partners would naturally be very short, and therefore these heavier 
odd multiplets would not contribute to the LOP relic density in any relevant way. 

Finally, in order to estimate the dark matter density, we shall restrict ourselves to the first level 
of odd KK modes, since they give the dominant contribution to the annihilation cross section. This 
is due not only to their relatively small masses with respect to the heavier modes in the tower, but 
also due to their larger couplings to the LOP. We have checked that the inclusion of the second 
KK level leads to a very small modification of the annihilation cross section and therefore of the 
freeze-out temperature and the predicted relic density. Let us finally comment that in the region 
consistent with the proper dark matter density, the mixing between the singlet and doublet particles 
is small and these particles lead to only a small contribution to the Higgs effective potential and 
the precision electroweak observables. 

Similar to the standard model fields discussed in the previous section, the boundary conditions, 
Eq. ( 1541) . lead to a set of equations which determine the masses of our odd multiplet fermions. For 
the KK modes that couple to the Higgs boson, in the case of vanishing Majorana masses for the 
odd fields, we find the following condition 



sm 



Xh 



f; 



h 



2 

+ S M2 S-m 2 = 0. (55) 



The solutions for the LOP are plotted in Fig. [3j The behavior of the LOP mass may be 
understood from the h = limit. In this limit, the singlet and doublet states don't mix and the 
singlet becomes the lightest odd particle for C2 < 0, while the doublet becomes the LOP for C2 > 0. 
At C2 = and h = 0, the singlet and doublet masses are degenerate. When h ^ the singlet and 
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Figure 3: The Dirac mass of the LOP, mi as a function of eg, the localization parameter for the odd fermions for 
three values of k = 1.5, 2.2 and 3.8 TeV. 

doublet states mix and their masses are split by the Higgs v.e.v. Only the lightest state mass is 
plotted in Fig. [31 In the presence of the Higgs, at C2 = 0, the LOP is an equal admixture of the 
singlet and doublet state. As we move away from c<i = 0, the roots of the determinant will start 
splitting into two clearly spaced masses. For c 2 < 0, the lighter mass is mostly a singlet state and 
the heavier one is mostly a doublet state. Since these are Dirac particles, the positive and negative 
roots of the determinant are equal. 

4.1 Odd Majorana Masses 

In the above, we have not considered the impact of Majorana mass terms that could in principle be 
written for this multiplet, both on the IR and the UV branes, as was done for the even neutrinos, 
and would modify the couplings to the Higgs boson. Including the Majorana masses for the odd 
multiplet, the equation determining the odd lepton masses is given by 



sin 



Xh 
fh 



1 2 

= 



S-Mz (Sm 2 ~ M IRo Sm 2 — e 2c2kL M UVo (^S^m 2 — M IRo S_ M2 

(56) 

The effect of introducing the Majorana mass terms can be seen in the different negative and 
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positive masses as one moves away from M[/y o ,M/R o = 0. Due to the different behavior of the 

functions S and S, by inspection, we can see that the positive and negative roots of Eq. (156|) will no 
longer be equal. The two Dirac states have been split into four Majorana states. These states can 
still be recognized as two mostly singlet and two mostly doublet states by comparing their masses 
to the charged states. In the following, we will sometimes refer to these states as singlet or doublet, 
where it should be understood that these are not really the original states, but those mixed by the 
Higgs. Without the mixing, the coupling between the singlet- singlet and the doublet-doublet states 
and the Higgs would vanish. Therefore, we expect the singlet-singlet coupling to be suppressed 
compared with the coupling of the mostly singlet state to the mostly doublet states. Additionally, 
looking at Fig. H] we see that as expected, the Majorana masses don't effect the mass of the LOP 
when the LOP is mainly a doublet state (purple circles for > 0), as the convergence of the 
different curves corresponding to the different values of Mm o and Muv a clearly show. 

The behavior of the LOP mass as well as its coupling to the Higgs may be studied by looking 
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at the roots of Eq. (1561) . Using the small z expansion for Sm, we obtain 



2 k I \ + c 2 



M IRo M UVo - e 2c2fci cos 



fh 



n 2 



2 -y e 2 C2fc L M/Ko + 4 (|±g M[/yo 



(57) 



which is valid for the case in which at least one of the Major ana masses is non- vanishing and c 2 < 0, 
in which the singlet becomes the LOP. This clearly shows that when only one of the Majorana 
masses is non- zero, we get a See-Saw effect governing the LOP mass. This behavior is confirmed in 
our numerical analysis plotted in Fig. HI 

Observe that if only the ultraviolet mass is non- vanishing, the mass of the LOP (which is mainly 
a left-handed singlet) is exponentially suppressed, unless c 2 — 0. Indeed, in the limit of vanishing 
infrared Majorana masses, Eq. (1571) reduces to 



z ~ — k - — c 2 



I \ e 2c2kL cos 



Xh 



2 7 2M, 



uv 



(58) 



As can be seen from Eq. ( 1581) . the higher operator coupling of the LOP to the Higgs is also suppressed 
for c 2 < 0. We will show that the annihilation cross section for a mainly singlet state is sufficiently 
enhanced only when the s-channel Higgs diagram becomes sizable and therefore, unless c 2 ~ 0, the 
Dark Matter density becomes very large compared to the experimentally observed value. 

As both Mj Ro and Mjjy o are turned on, we see an abrupt change in the behavior of the mass 
spectrum for c 2 < 0, which becomes independent of the exact value of M UVo and only depends on 
k and M IRo , 

k M IRo 



z ~ 



2-°*) ■ ( 59 ) 



The LOP mass in this case is of the order of k, and does not show an explicit dependence on the 
Higgs vaccum expectation value. Indeed, as can be seen from Eq. ( IBTI) . the effective coupling to the 
Higgs for c 2 < is exponentially suppressed. Therefore, as happens in the case of vanishing Mi Ro , 
a good dark matter candidate may only be obtained for values of c 2 > 0. 
When c 2 ~ 0, the mass can be approximated by: 



„M IRo M UVo -e 2c * kL cos 
z ~ k- 



Xh 

.fh 



e 2c * kL M IRo + 4M ; 



UVo 



(60) 



We see that for Mi Ro , M\jy o ~ 0(1), for values of h in the linear regime and very small values of 
c 2 , we can get a cancelation resulting in very small LOP masses. This behavior is clearly portrayed 
in Fig. HI 

Finally, let us analyze the case M IRo ^ and M UVo = 0. As M IRo is turned on, it strongly 
modifies the spectrum with respect to the Dirac case for negative values of c 2 . In this case, the 
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LOP becomes mostly a right-handed singlet and its mass is given by 



cos ¥ 

z „ -fc(l+2c 2 ) y hi . (61) 
M IR 

As seen in Eq. ( IBTj) . the LOP mass in this case is, again, of the order of the weak scale but with an 
explicit Higgs v.e.v. dependence induced by a higher order operator coupling with a characteristic 
scale of the order of the KK masses. Therefore the coupling of the LOP to the Higgs becomes 
sizable for KK masses of the order of the TeV scale, allowing for the possibility of a dark matter 
candidate for ci < 0. 

In the following section, we shall perform a more precise quantitative analysis of the masses and 
couplings associated with the annihilation cross section of the singlet state for both the Majorana 
and Dirac cases. 



4.2 Couplings of the Odd Leptons 
4.2.1 Higgs Couplings 

To calculate the couplings of the Majorana and Dirac states with the Higgs, the profile function 
of the odd leptons which couple to the Higgs bosons need to be computed. The mass eigenstate 
profile functions are given in terms of combinations of the profile functions without the Higgs. In 
the particular case of the neutral odd leptons these are admixtures of the neutral states belonging 
to the bidoublet and the singlet state, with normalization coefficients C2, C3 and C5. 

The boundary conditions determine the coefficients C3 and C5 as functions of CV The fermion 
profile functions in the presence of the Higgs are given by: 
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(62) 
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Figure 5: Higgs couplings to the Dirac particles for three values of k ~ 1.5, 2.2 and 3.8 TeV corresponding to 
mi ~ 2, 3 and 6 TeV for C2 = from left to right, as a function of the singlet LOP mass mi. 
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where the functions S and 5 are functions of x§ and the masses z of the odd fermions. 

For the doublet and singlet states mixed by the Higgs, a non-trivial solution may be only obtained 
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Figure 6: Higgs couplings to the Majorana particles for three values of k ~ 1.5, 2.2 and 3.8 TeV corresponding to 
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if the following relations are fulfilled. 
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For the neutral leptons which decouple from the Higgs, instead, the following relations are fulfilled: 

C 3 = -C 2 (70) 







(71) 



This implies that only the symmetric combination of neutral bidoublet states with coefficients 
given by Eq. (|68l) and (|69l) . couple to the Higgs. Moreover, the normalization coefficient Ci may be 
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Figure 7: Higgs couplings to the Majorana particles for three values of k ~ 1.5, 2.2 and 3.8 TeV as a function of 
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computed by demanding well normalized functions, namely, 

C 2 = Ei ^ { ^ + lfkl2) dx) V2 . (72) 

The above definition is appropriate in the Majorana case, in which the left-handed components of 
the fermions acquire contributions from both the original left-handed and the (charge conjugate) 
right-handed modes, Eqs. (!62j) - (l67|) . In the Dirac case, the left-handed and right-handed functions 
acquire equal normalizations and therefore the proper factor C 2 is equal to the one computed above 
divided by \/2. In the following, we will keep the above definition of C 2 for both the Majorana and 
Dirac cases and take care of the proper a/2 factors explicitly 

To calculate the couplings, we also need the Higgs profile and normalization: 

fh = C h e 2kx (73) 

C h = 95 (74) 

v/o a ( x ) 2dx 
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Defining H(mj, rrij): 



Eim^mj) = - e —C h C* 2 (m t )C 2 ( mj )f h [f R * (m.) [fl( mj ) + f L ( mj )] - [fifa) + /l*(m*)] /*(m,)] 

(75) 

the left-right couplings of the Higgs with the different states, ty l L Htyi , in the Majorana and Dirac 
cases can be written as: 

Xfj = / (E(mi,mj) + E*(mj,mi)) dx (76) 
Jo 

Ag = 2 f E{m h rrij) dx (77) 
Jo 

the factor of 2 in the Dirac coupling is due to the definition of the C 2 factor discussed above. 
Observe that, while in the Majorana case Xf^j = A^f *, there is no such relation in the Dirac case. 

The different couplings are plotted in Figs. El -[3 For the Dirac case and c 2 < 0, represented in 
Fig. small values of the masses are obtained for smaller values of c 2 . The left-right couplings of 
the singlet and doublet states, Af 2 acquire large values for negative values of c 2 . This stems from 
the fact that for this case, the left-handed singlet component is localized towards the IR brane. 
As c 2 goes to 0, the localization effects become less pronounced and this coupling starts getting 
suppressed. The Af^ couplings have the opposite behavior to the cross couplings. For c 2 negative, 
the mass difference between the singlet and the doublet state is large, while their mixing is small. 
Since the self-couplings of the mass eigenstates are induced by the product of the singlet and doublet 
components of these states, they become very suppressed. However, as c 2 goes to zero, the mass 
eigenvalues become symmetric and antisymmetric combinations of the singlet and doublet states, 
and the self couplings of the mass eigenstates become large, while the cross couplings tend to zero. 

In the Majorana case with M\jy o = 0, although the quantitative values are not the same, the 
Afi and A^ 4 couplings behave similarly to the Af^ couplings, since mi and are the two mostly 
singlet states, which are split due to the non-zero Mir o . When both the Majorana masses are non- 
zero, we see an abrupt change in the behavior of the couplings. As mentioned before, for c 2 < the 
self-coupling of the lightest state is exponentially suppressed. This behavior is clearly demonstrated 
in the Higgs dependent part of the approximation for the LOP given in Eq. (1571) . The couplings of 
the LOP to the mostly doublet states, however, continue to be large. 

4.2.2 Couplings to the Z and W ± Bosons 

The Z couplings to the lepton sector are defined in a similar manner to the couplings with the 
quark sector [25], [26]. However, in the lepton case Qx = and the neutral states that couple to 
the Higgs have C 3 = C 2 , implying that f1 R {h) = fi R {h). Therefore, the Z, and all its KK modes, 
as well as the neutral components of the SU(2) R gauge bosons don't have any couplings with any 
two of these states. However, the orthogonal neutral state in the bidoublet, which does not couple 
to the Higgs has C3 = — C 2 and C5 = 0. Hence, an off-diagonal coupling exists between this mode, 
the neutral states that couple to the Higgs and the Z. 
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Figure 8: M /± couplings to the Dirac particles for c-i < and three values of k ~ 1.5, 2.2 and 3.8 TeV, which, for 
C2 = 0, correspond approximately to ui\ ~ 2, 3 and 5.6 TeV, from left to right as a function of the singlet LOP 
mass. Larger values of mi are associated with larger values of c%. 

The W ± couple the charged fermions with the neutral components. In component form, the 
coupling is between f£ji{h) and (h). The profile functions and their normalization coefficients 
for the neutral components were given in Eqs. flBT|) - (|68j) . The charged fermions and the neutral 
component of the bidoublet which don't couple to the Higgs state are both governed by the same 
five dimensional wave-function, namely: 

f l L = C i e^ 1+ ^ kx L Ma (78) 
f R = C^^S.m, (79) 

These fermion masses are given by the roots of S-m 2 an d since the Majorana masses don't influence 
them, they are always Dirac states. Further, it can be shown trivially that the W~ coupling to the 
LOP and the positive charged state is equal to the coupling of the LOP to W + and the negative 
charged states. In the annihilation cross-section, we will only be interested in the couplings between 
the two charged fermions and the LOP. We will denote these couplings by gw L R - The expression 
for these couplings is given in Appendix B. Similarly, for the Z we will only be interested in the 
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Figure 9: W ± couplings of the Majorana particles for C2 < and three values of k ~ 1.5, 2.2 and 3.8 TeV, which, 
for C2 = 0. correspond to mi ~ 1.7, 2.5 and 4.7 TeV from left to right as a function of the singlet LOP mass. Larger 
values of mi are associated with larger values of ci. 

couplings between the LOP and the N' state, the neutral bidoublet component that does not mix 
with the Higgs, and we will denote these couplings by gz L R ■ 

The gauge boson couplings are plotted in Figs. M - EH Again we see that the Dirac and the 
Majorana, M\jy o = couplings behave in a similar way. The coupling of the mostly singlet state 
to the charged or neutral fermion is obtained through the mixing with the bidoublet states. As 
discussed before, this mixing is small, for C2 < 0, and increases for larger values of 0%. As C2 
approaches 0, in the Dirac case, we expect that since the mixing is maximal, the couplings 
should approach the neutrino- lepton SM coupling, g w /\^2, reduced by a factor l/\/2> due to the 
mixing of the singlet with the doublet state, times another factor l/v2 due to the projection of the 
neutral doublet state on the SU(2) partner of each of the charged fields. This behavior is clearly 
seen in Fig. El where only values of C2 < are plotted and increasing values of mi are associated with 
larger values of C2- For C2 < 0, the left-handed states which are located towards the infrared brane 
couple more strongly to the charged states than the right-handed states. Moreover, the couplings of 
the Z and the W ± become proportional to each other, with a coefficient of proportionality governed 
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Figure 10: couplings to the Majorana particles for three values of k ~ 1.5, 2.2 and 3.8 TeV as a function of 
C2 for the singlet LOP. The smaller values of k correspond to the smaller couplings (the lower curve). 



by cos 6^, namely 

9z L , R = ±~ (80) 

cos Vw 

The additional factor of y/2 that appears between SM couplings of neutrinos to the Z and W, is 
not present in this case. 

In the case of zero ultraviolet Majorana mass but non- vanishing Mm , depicted in Fig. [HI 
the behavior is similar to the Dirac case but the couplings are reduced due to the larger singlet 
components of the Majorana particles. Also, there is a sizable reduction of the left-handed couplings 
due to the larger right-handed component of the Majorana state. Observe that as M\jy is turned 
on, for c 2 < 0, the couplings to the gauge bosons become independent of c 2 . 

Finally, as c 2 becomes positive, the couplings increase, due to the larger bidoublet component 
of the LOP. 
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Figure 11: Feynman Diagram for the process N i + Ni — + t + t 



4.3 Annihilation Cross Section 

We will denote the neutral states mixed by the Higgs by Ni, where % = 1,2, or i — 1,2,3,4 for 
the two Dirac or four Majorana states respectively, where i labels the states in increasing order of 
their absolute masses. C ± will denote the charged fermions and, as said before, N' will denote the 
bidoublet neutral fermion which does not couple to the Higgs. The Ni is the LOP, our dark matter 
candidate. 

Ignoring co- annihilation effects, we consider the following five dominant processes for N\Ni 
annihilation: Ni + Ni — - > t i , H H, Z Z and W + W~ (observe that due to the cancelation of the 
Z coupling to the states that couple to the Higgs, the Z H annihilation channel is suppressed). 
The Feynman diagrams contributing to each of these processes are shown in Figs. [TT] - [TH The 
virtual iVj exchanges in these diagrams should be understood to be summed over i, where i as noted 
above runs over the appropriate index depending on whether we are considering the Dirac or the 
Majorana case. The v in the following formulae is the relative velocity between the initial particles 
in the center of mass frame. Xmt, ^hzz, ^hww an d Xh are the couplings of the Higgs to the top, 
the W ± and Z bosons, and itself, which were discussed in Refs. (23] and [25] . 

4.3.1 iVx + iVx t + t 

Due to the cancelation of the coupling of N% to the Z, the annihilation into fermion pairs proceeds 
via an s-channel Higgs interchange, and is therefore proportional to the corresponding fermion mass. 
Therefore, only the top contributes in a relevant way. The Dirac and Majorana cross-sections are 
given by the same formula, but the Higgs coupling should be understood to be the one appropriate 
for each case. Assuming mi > m t , we obtain: 




(81) 
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Figure 12: Feynman Diagrams contributing to the process Ni + N\ — > H + H 
4.3.2 N x + ^ -> if + H 

The annihilation into Higgs pairs proceeds via an s-channel Higgs interchange diagram, which is 
subdominant, and the t-channel interchange of the neutral odd fermions. The result, in the limit 
ran < Cmi,m 2 , is given by: 



< ™ >HH = 




where we have taken the couplings to be real. The sum runs over the two Dirac states labeled by 
their masses, mi and m 2 . In the case of real couplings, the Majorana cross-section can be simply 
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seen from the above with the replacement Xf^ = 1/2 (Afj + Aj^), and the indices now run over the 
four Majorana states: 
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(83) 



4.3.3 Ni + Ni^W + W~, Z + Z 

The annihilation into the and Z gauge bosons also proceeds via the s-channel interchange of 
a Higgs, plus the t-channel interchange of the charged fermion and the neutral fermion N', 
respectively. In the formula below, the label G corresponds to either the W ± or Z gauge bosons, 
and a = 1 for the W + W~ cross-section and a = 1/2 for the ZZ case. The diagrams contributing to 
the process in the Dirac case are given in Fig. [131 For the Majorana case, two additional diagrams 
contribute, and are given in Fig. [TH Using the properties of the Majorana couplings, one can 
demonstrate that these new diagrams are equal to the ones associated to the cross diagrams in 
the amplitudes for the annihilation into the W ± (due to the interchange of the fermion of opposite 
charge) and Z gauge bosons. Therefore, the cross-section is given by the same formula for both the 
Dirac and Majorana cases, but with appropriate couplings, and a factor (3 = 2 for the Majorana 
case, and (3 = 1 for the Dirac case. Although in the numerical analysis the full annihilation cross 
section was used, for simplicity, we will only quote the cross-section for the longitudinal modes, in 
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Figure 13: Feynman Diagrams contributing to the process Ni + N% — > W + W ,Z + Z. The intermediate state is 
either the charged fermion C for the case, or the orthogonal bidoublet, N' for the Z. 
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These cross-sections are plotted in Figs. [ToTfTTl For negative Ci (smaller values of mi), we observe 
an interesting correlation between the annihilation cross sections into W^, Z and Higgs pairs. We 
see a dominance of the longitudinal modes for this range of values of C2 and the magnitudes of 
the W ± ) Z Z and H H cross-sections obey the 2:1:1 behavior expected due to the Goldstone 
equivalence theorem. For larger values of C2, the bidoublet component of the LOP increases and 
the transverse components of the gauge bosons are no longer subdominant in their contribution to 
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Figure 14: Additional Feynman Diagrams contributing to the process N± + N± — > W + W , Z + Z for the Majorana 
case. The intermediate state is either the charged fermion C for the W case, or the orthogonal bidoublet, N' for 
the Z. 

the annihilation cross section. 

Our extensive numerical and analytic study showed that for negative C2, the major contributions 
to the W and Z Z cross-sections are due to the s-channel Higgs exchange. In the H H cross-section, 
this is matched by the contribution from the virtual exchange of the iVj in the t-channel. Therefore, 
as emphasized before, for ci < 0, a sizable annihilation cross section may only be obtained when 
the Higgs coupling to the LOP becomes of order one. 

4.4 Dark Matter Density 

We shall follow the standard formalism for the calculation of the thermal dark matter density [44], [45] 
In calculating the annihilation cross-sections, we used the non-relativistic approximation for the ini- 
tial particles. The cross-section used in calculating the dark matter density is the sum of all the 
different contributions in the previous sections and will be denoted by < av >t, and xp = m/T F 
as usual. The relative velocity is related to the freeze-out temperature: 



The non-relativistic expansion of the thermal annihilation cross section may be expressed as 
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(85) 
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< av >t= Co + °\ < v 2 >— °o + 6 (Ti/xp. 



(87) 



The dark matter density is then given by 




(88) 
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Figure 15: The cross-section contributions to the annihilation of the Dirac LOP for (from top to bottom) k — 1.5, 
2.2 and 3.8 TeV. 

where c = 1/2, M Pl = 1.2 x 10 19 GeV, g* = 112, s = 2889.2/cm 3 , p c = 5.3 x lCr 6 GeV/cm 3 , g = 2 
is the degrees of freedom of our dark matter candidate and 7 = 2 or 1 to account for the antiparticles 
for the Dirac and Majorana case respectively. We take g$ = 2 in the Dirac case since we compute 
the density of the particle and antiparticle separately. The factor 7 = 2 in the relic density then 
accounts for the duplication of the density from both the Ni particle and the antiparticle in this 
case. 

To check the veracity of our calculations, we extensively studied the limit of the Majorana case 
reducing to the Dirac case as M IRo and M UVo go to 0. As the Majorana masses go to 0, the two singlet 
masses start getting degenerate in mass, and the N\ N\ and the N2 cross-sections become equal. 
To properly analyze this limit, then, we must take co-annihilation between the lightest Majorana 
sates into account [45J. One can check that the coupling of the Higgs to each of the degenerate 
LOP Majorana states, HNiNi, becomes equal to the one of the Higgs to the LOP in the Dirac case. 
Moreover, the cross coupling of the Higgs, HN1N2, vanishes identically in the limit of vanishing 
Majorana masses. Further simple relations exist between the Higgs couplings to fermions in the 
Majorana and Dirac cases. One can check that due to these relations the annihilation cross section 
into Higgs states in the Majorana case NiNi — > H H become the same as the NiN\ — > H H cross 
section in the Dirac case. The same happens in the case of annihilation into fermions. 

In the case of gauge bosons the situation is more complicated. As noted in calculating the 
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Figure 16: The cross-section contributions to the annihilation of the Majorana LOP with Mm o = 0.5 and Mjjv = 
0, for (from top to bottom) k = 1.5, 2.2 and 3.8 TeV. 

couplings, the gauge boson couplings in the Majorana case are reduced by a factor l/y/2. This 
implies that for the W ± and Z Z, the interference between the t and s-channel diagrams is reduced 
by 1/2 and the t-channel diagrams by 1/4 compared to the Dirac cross-section case. These factors 
are exactly compensated for by the extra diagrams that contribute in the Majorana case (Fig. [T4"]) . 
and therefore one obtains that the equality of annihilation cross sections defined above for the Higgs 
final states extends to all final states. We also verified that the N\ N 2 annihilation cross-section is 
exactly in this limit. Including co-annihilation between the two Majorana states |45j . the effective 
degrees of freedom are then 4, and the effective cross-section is | a D , where a D is the annihilation 
cross section between Ni and its antiparticle in the Dirac case. Therefore, in this limit the dark 
matter density due to the Dirac and Majorana cases in the limit Muv ,ir ^ is exactly the same, 
as expected. 

We plotted the values of c 2 and mi leading to the correct dark matter density in both the Dirac 
and the Majorana cases in Fig. [18j We restricted the values of k > 1.2 TeV, for which the SM Higgs 
couplings are close to their SM values. This provides a lower bound on the LOP mass and on the 
value of C2- The bands in this figure represent the relic density uncertainty. We see that in the 
Dirac case, we can only obtain the correct dark matter density for values of 2 TeV £ mi > 1 TeV 
and C2 > —0.4, consistent with the exchange symmetry. The value of k is correlated with the value 
of mi, and is constrained to be in the range 2 TeV Z k > 1.2 TeV for the above range of masses. 
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Figure 17: The cross-section contributions to the annihilation of the Majorana LOP with Mm o = 0.5 and Mjjv = 
0.5, for (from top to bottom) k = 1.5, 2.2 and 3.8 TeV. 



In the Majorana case, for M.jjy = 0, we are able to obtain the correct dark matter density 
for a larger range of mi values, which extend from about 3 TeV up to the lowest values of about 
mi ~ 500 GeV, and for the range of negative values of c 2 > —0.4. The values of mi are, again, 
correlated with the values of k, which is in the range 3 TeV £ k £ 1.2 TeV. For Muv 0> instead, 
we see that we can only obtain the correct dark matter density for C2 ^ 0. Due to the effect of the 
Majorana masses, the singlet LOP state mass is still significantly smaller than the doublet mass in 
this case and therefore co-annihilation effects may be ignored. 

The values of C2 obtained for the case of vanishing ultraviolet Majorana masses are fully compat- 
ible with the identification of the LOP with the odd partner of one of the right-handed neutrinos. 
Indeed, as can be seen from Fig. [21 for values of C\ ~ 0.5-0.7, the proper neutrino masses are ob- 
tained for values of —0.4 < c 2 < —0.1, for which a proper dark matter candidate may be obtained, 
with a mass 0.5 TeV < mi < 2.5 TeV. 

The situation is different for non- vanishing values of the ultraviolet Majorana mass Muy o . In this 
case, a proper dark matter candidate may only be obtained for values of c 2 ~ 0. Such values of c 2 
are incompatible with the exchange symmetry if the ci's of the three generations are approximately 
the same, as assumed in this article. Therefore, a proper dark matter candidate would demand 
that the Majorana mass at the UV brane is either zero or smaller than exp(2c2fc L), since in such 
a case, according to Eq. (1571) . masses of the order of the weak scale and couplings to the Higgs of 
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Figure 18: Parametric plot of mi, the mass of the LOP, versus C2, the localization parameter of the odd fermions, 
when Qdm ~ 0.23 ±0.1. The two lines for each value of the Majorana masses are associated with the upper and 
lower bound on Qdm- 

order one would be obtained. 

Observe that the difficulty in obtaining a proper dark matter candidate for Muy o ^ stems from 
the fact that we have assumed equal bulk mass parameters for the fermions interchanged by the 
Z 2 exchanged symmetry introduced in Refs. [13]. Although this is an attractive possibility, since it 
allows a connection between the dark matter properties and the neutrino masses, this does not need 
to be the case. Values of c 2 > would be allowed in the more general case, or for a more general 
discrete symmetry. Alternatively, if the value of c\ of the left-handed leptons was allowed to be 
different for the three generations, values of c 2 > would be consistent with values of c% compatible 
with the relatively small electron mass. 

4.5 Direct Dark Matter Detection through Higgs Exchange 

The annihilation cross section for the odd neutrinos becomes of the proper size for sizable values of 
the coupling of the odd neutrino to the Higgs boson, An ~ 0.3 - 0.7, with larger couplings associated 
with larger LOP masses. Such large couplings induce a relatively large scattering cross section of 
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the odd lepton with nuclei that may be probed at direct dark matter detection experiments. More 
quantitatively, the spin-independent elastic scattering cross-section for an odd lepton scattering off 
a heavy nucleus is: 



4m 
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a SI = ^(Zf p + (A-Z)f n ) 2 (89) 
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where m r = mN+m ^ and is the mass of the nucleus. The factors / PjTl are given by 
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where the quark form factors are = 0.020±0.004, = 0.026±0.005, = 0.118±0.062, f TG « 
0.84, = 0.014 ± 0.003, = 0.036 ± 0.008, /£ = 0.118 ± 0.062 and 0.83 [46J. Hence, we 

find that the contribution is 
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where we have neglected the differences between the proton and the neutron mass and have used 
the fact that the neutron and proton fa factors are relatively similar. Assuming that the mass of 
the odd neutrino is much larger than that of the nucleus we have m r ~ m N ~ Am p and 
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where <7s//v4 4 is the neutrino-nucleon spin-independent cross-section. From Eq. (195j) . the spin- 
independent cross-section scales as A^/mf^ and therefore direct dark matter detection experiments 
like CDMS can put strong constraints on regions of small mn and large An. 

As discussed above, in the Majorana case, an odd neutrino with a mass of about 700 GeV and a 
coupling to the Higgs of about 0.35 leads to an acceptable dark matter density. The spin independent 
cross section obtained in such a case for a Higgs mass, m# ~ 130 GeV, is about 1.4 x 10~ 43 cm 2 . 
The current limit coming from CDMS, from a combination of the Ge data and under the standard 
assumptions of local dark matter density distribution, is about 2.5 xlO -43 cm 2 [17]. The XENON 
experiment puts a slightly weaker limit for this range of LOP masses |48] . Therefore, the predicted 
spin independent cross section is only a factor of a few lower than the current experimental limits. 
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For larger masses, of about 1 TeV, the Higgs couplings grows to about 0.38 and the predicted cross 
section is therefore about 1.65 x 10~ 43 cm 2 , while the CDMS bound is about two times larger. In 
the Dirac case, the couplings are about fifty percent larger than in the Majorana case and therefore 
the predicted cross section for a mass mi ~ 1 TeV is slightly above the CDMS reported bound. One 
would be able to conclude that the Dirac case for ma ~ 130 GeV is therefore disfavored. There are, 
however, uncertainties of order of a few associated with the local dark matter density distribution 
and the nuclear form factors which should be taken into account before ruling out a specific model. 
It is expected that both the XENON and CDMS experiments will further improve their sensitivity 
by about an order of magnitude by the end of 2009 [39], [50]. Therefore, even considering possible 
uncertainties associated with the local density and the nuclear form factors, the minimal model 
discussed in this article should be probed by these experiments in the near future. 

Let us comment that the mass of the N\ particle may be in the appropriate range to provide 
an explanation of the anomalous excess in electrons and positrons observed by the Pamela [51] 
and ATIC [52] experiments. However, since in these model these particles decay mostly into Higgs 
and gauge bosons, if these particles are distributed throughout the halo of the galaxy, an excess of 
positrons of the size observed by these experiments will need a large boost factor enhancement and 
would probably lead to an unobserved excess of antiprotons [53j,[54j. 

5 Conclusions 

In this article, we have considered the question of incorporating the charged and neutral leptons 
in a Gauge-Higgs Unification scenario based on the gauge group SO(5) x U{l)x in warped extra 
dimensions. These models are attractive since the SO{4) = SU{2)l x SU{2)r subgroup of SO(5) 
incorporates in a natural way the weak gauge group as well as an appropriate custodial symmetry 
group in order to suppress large contributions to the precision electroweak observables. Moreover, 
the fifth dimensional components of 5*0(5)/ SO (4) gauge bosons have the right quantum numbers 
to play the role of the Higgs doublet responsible for the breakdown of the electroweak symmetry. 

We have shown that, similar to the quark sector, the leptons can be incorporated by including 
the left-handed zero modes in a fundamental representation of SO(5) and the right-handed charged 
leptons in a 10 of SO (5). The model includes right-handed neutrinos which are singlets under the 
SO (4) x U(l)x group and can therefore acquire localized Majorana masses on the IR and UV branes. 
The simple inclusion of the right-handed neutrinos in the same multiplet as the charged leptons, 
fails to produce the correct lepton masses. The correct charged lepton and neutrino masses may 
be obtained from a three multiplet structure similar to the quark case. The bulk mass parameters 
ci, C2 and C3 of the left-handed leptons, right-handed neutrinos and right-handed charged leptons, 
respectively acquire values c\ ~ 0.5 - 0.7, —0.4 < c% £ —0.1 and —0.9 < C3 < —0.5, where larger 
negative values of C3 correspond to the first generation leptons. 

We have further investigated the possibility of incorporating a dark-matter candidate by in- 
cluding an exchange symmetry, under which all SM leptons multiplets are even, and which ensures 
the stability of the lightest odd lepton partner. We therefore analyzed the possibility of associat- 
ing the dark matter with the lightest neutral components of the odd leptons, transforming in the 
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fundamental representation of SO (5). We have shown that these neutral components have inter- 
esting properties. The neutral leptons that couple to the Higgs do not have self couplings to the 
Z-boson. However, these neutral states couple to the orthogonal combination of neutral states in 
the bidoublet and to the Z, as well as to the charged leptons and the W-gauge boson. 

We computed the couplings of the neutral odd leptons to the gauge bosons and the Higgs in 
a functional way and computed the dominant contributions to the annihilation cross section into 
Higgs, neutral and charged gauge bosons and fermions (top-quark) final states. We considered the 
cases in which the Majorana masses for the neutral leptons vanish in both branes (Dirac case) as 
well as the case in which at least one of them is non- vanishing. By doing that, we have shown that 
in the Dirac case, a proper dark matter candidate may be obtained for masses of about 1 TeV to 
2 TeV and localization parameter —0.4 < c 2 < —0.1 in agreement with the exchange symmetry. 
If only the Majorana mass in the infrared brane is non-vanishing, one obtains lower values of the 
required odd-lepton masses, which may be of about 500 GeV to 2.5 TeV, and a range for the 
localization parameter —0.4 < c 2 < —0.1. Finally, in the case that the Ultraviolet Majorana mass 
is different from zero, the self-couplings of the LOP with the Higgs is exponentially suppressed for 
c 2 < and becomes non- vanishing only for c 2 — 0, for which a proper dark matter may be obtained 
with a mass similar to the one obtained when only Mm is different from zero. This last possibility 
is incompatible with the exchange symmetry and the proper neutrino masses if a common values 
of ci is assumed for the three families. 

The collider signatures of these models have been previously discussed in the literature. The 
odd leptons introduced in this article will be hard to detect at collider experiments, since the masses 
of the charged and neutral non-LOP odd leptons are above a few TeV, and they have relatively 
weak interactions. In all cases, a proper dark matter candidate is obtained for values of the self- 
coupling of Ni to the Higgs of about 0.3 - 0.7, with larger Higgs couplings corresponding to larger 
LOP masses, for which the cross section of the dark matter with nuclei becomes sizable. We have 
computed the dark matter cross section with nuclei and have shown that these models will be 
probed by the CDMS and XENON direct dark matter detection experiments in the near future. 
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APPENDIX 



A Profile Functions at h = 0. 

In the h = gauge, we redefine ip = a 2 {x§)il) and we write our vector-fermionic fields in terms of 
chiral fields. We can KK decompose the fermionic chiral components as, 



L.R 



n=0 



where / is normalized by: 



dx 5 a 1 {x b )f n f m = 5 n 



(A.l) 



(A.2) 



Therefore the profile function for the zero mode fermion corresponds to a~ 1 l 2 (xs)fQ. 

From the 5D action, concentrating on the free fermionic fields, we can derive the following first 
order coupled equations of motion for fL,R,n, 



(d 5 + M)f R>n = (z/a(x 5 ))f L ,n, {d 5 - M)f L ,n = -{z/a{x 5 ))l 



R 



(A.3) 



We see from Eq. ( 1A.3I) that we can redefine Jr } l,', 



-Mx f 



fR,L,n an d relate the opposite chiral 



component of the same vector-like field by f^ n = (—a(x 5 )/z)d 5 fL, n - For the left handed field having 
Dirichlet boundary conditions on the UV brane, we can derive a second order equation for the chiral 
component f^ n : 



2M d 5 + 



f L ,n = 



(A.4) 



the solution of which we shall call Sm(x 5 , z), with boundary conditions Sm (0, z) = 0, S' M (0, z) = z. 
Similarly, if the right-handed field fulfills Dirichlet boundary conditions on the UV-brane, we 



can redefine fR^,n = e s fR,L,n an d then relate the opposite chirality via f^, 
We can further write the equation of motion for /p n : 



{a(x 5 )/z)d 5 f_ 



R,n- 



2M)d 5 + 



0. 



(A.5) 



We shall correspondingly denote the solution to this equation with 5*_m(^5, z), fulfilling the bound- 
ary conditions S-.m(Q,z) = 0, S'_ M (0,z) = z. 
The solution to Eq. (IA.4I) is given by [3]: 



at \ 7TZ -c- 



^5 



+c 



1 1 Y, 



+c 



ka(x$) 



Yi 



+c 



ka(x 5 ] 



(A.6) 



The solution for Eq. (1A.5j) . S-m, is given by Eq. flA.6j) with the replacement c 



-c. 
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B Coupling of the charged gauge bosons 

Following the notation of Ref. |26j, the boson profile functions are given by: 



fUh) 





'Xh 


(cos 







Of + -= si 
1 V2 



sm 



Xh 

fh . 



( ) --j=C(x s ) sin 



Xh 

fh 



(B.l) 



fc L (h) 



- (s{x s ) ( (cos W - l) C? R + V2 sin [^] Cf ) + CM (l + cos [^1 ) cA 



(B.2) 



S(x 5 ) cos 



+ C(a; 5 ) 1 - cos 



Xh 

fh 



(B.3) 



The normalization coefficients , in terms of C9 are given by: 

1,1k 1l 



C 



G 

L 



—4 cos 


Xh 
Jh_ 


C{L)'S{L)' + C h he 2kl 


' sin 


Xh 
_fh_ 


(C(L)'S{L) + C{L)S{Ly) 


y/2S(L)' (c h he 2kL cos 


Xh 
J". 


l S(L)- 


f 2 sin 


Xh' 

fh 





r c C(L)' 
^S{L)> 



(B.4) 
(B.5) 



The five dimensional weak coupling is defined as g 5w = g w V~L. In terms of these, the Dirac 
couplings for W + and Z, (denoted by G), are given by: 



-9B W J o (fl: R (h) (ft R (h)fh(h) + fUWaW - fl R (h)f G R (h))) dx 5 



(B.6) 



Again, due to our choice of normalization for the coefficient Cg, the factor y/2, coming from the 
definition of the W ± fields, is not present in this expression. The Majorana couplings are given in 
terms of the Dirac couplings, 

(B.7) 



with g D given in Eq. ( IB. 61) . 



42 



References 

[1] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) |arXiv:hep-ph/9905221| . 

[2] H. Davoudiasl, J. L . Hewett and T. G. Rizzo, Phys. Lett. B 473, 43 (2000) 
[arXiv:hep-ph/9911262| . 

[3] A. Pomarol, Phys. Lett. B 486, 153 (2000) |arXiv:hep-ph/9911294| . 

[4] S. Chang, J. Hisano, H. Nakano, N. Okada and M. Yamaguchi, Phys. Rev. D 62, 084025 (2000) 
[arXiv:hep-ph/991249~8j . 

[5] T. Gherghetta and A. Pomarol, Nucl. Phys. B 586, 141 (2000) [arXiv:hep-ph/0003129| ; 

[6] C. Csaki, J. Erlich and J. Terning, Phys. Rev. D 66, 064021 (2002) | arXiv:hep-ph/0203034] . 

[7] J. L. Hewett, F. J. Petriello and T. G. Rizzo, JHEP 0209, 030 (2002) |arXiv:hep-ph/0203091| . 

[8] H. Davoudiasl, J. L. Hewett and T. G. Rizzo, Phys. Rev. D 68, 045002 (2003) 
|arXiv:hep-ph/0212279~| . 

[9] M. S. Carena, E. Ponton , T. M. P. Tait and C. E. M. Wagner, Phys. Rev. D 67, 096006 (2003) 
[arXiv:hep-ph/0212307| . 

[10] M. S. Carena, A. Delgado, E. Ponton , T. M. P. Tait and C. E. M. Wagner, Phys. Rev. D 68, 
035010 (2003) larXiv:hep-ph/0305188] . 

[11] M. S. Carena, A. Delgado, E. Ponton , T. M. P. Tait and C. E. M. Wagner, Phys. Rev. D 71, 
015010 (2005) larXiv:hep-ph/0410344] . 

[12] K. Agashe, A. Delg ado, M. J. May and R. Sundrum, JHEP 0308, 050 (2003) 
|arXiv:hep-ph/03080~3~6~| . 

[13] K. Agashe, R. Cont ino, L. Da Rold and A. Pomarol, Phys. Lett. B 641, 62 (2006) 
[arXiv:hep-ph/060534T] . 

[14] M. Carena, E. Ponto n, J. Santiago and C. E. M. Wagner, Nucl. Phys. B 759, 202 (2006) 
|arXiv:hep-ph/0607i06| . 

[15] M. S. Carena, E. Pont on, J. Santiago and C. E. M. Wagner, Phys. Rev. D 76, 035006 (2007) 
[arXiv:hep-ph/0701055| . 

[16] N. Arkani-Hamed and M. Schmaltz, Phys. Rev. D 61, 033005 (2000) | arXiv:hep-ph/9903417| . 
[17] Y. Grossman and M. Neubert, Phys. Lett. B 474, 361 (2000) |arXiv:hep-ph/9912408| . 
[18] S. J. Huber and Q. Shan, Phys. Lett. B 498, 256 (2001) j|arXiv:hep-ph/0010195] . 



43 



[19] N. S. Manton, Nucl. Phys. B 158, 141 (1979); Y. Hosotani, Phys. Lett. B 126, 309 
(1983) ; H. Hatanaka , T. Inami and C. S. Lim, Mod. Phys. Lett. A 13, 2601 (1998) 
[arXiv:hep-th /98050 67|; I. An toniadis, K. Benakli and M. Quiros, New J. Phys. 3, 20 
(2001) |arXiv:hep-th/0108005|; M. Kubo, C. S. Lim and H. Yamashita, Mod. Phys. Lett. 
A 17, 2249 (2002) [arXiv:hep-ph/0111327|; C. Csaki, C. Grojean and H. Murayama, Phys. 



Rev. D 67, 085012 (2003) |arXiv:hep-ph/ 0210133|; N. Haba, M. Harada, Y. Hosotani 
and Y. Kawamura, N ucl. Phys. B 657, 169 (2003) [Erratum-ibid. B 669, 381 (2003)] 
|arXiv:hep-ph/0212035|; C. A. Scrucca, M. Serone and L. Silvestrini, Nucl. Phys. B 669, 128 
2003) [arXiv:hep-ph/030"4220] ; G. Burdman and Y. Nomura, Nucl. Phys. B 656, 3 (2003) 
|arXiv:hep-ph/02102"57[ ~ 

[20] R. Contino, Y. No mura and A. Pomarol, Nucl. Phys. B 671, 148 (2003) 
[arXiv:hep-ph/03062'59| . 



[21 

[22; 

t 23 
[24 

[25; 
[26 
[27 

[29 
[30 
[31 
[32; 
[33 
[34 



K. Agashe and R. Contino, Nucl. Phys. B 742, 59 (2006) [arXiv:hep-ph/0510164| . 

R. Contino, L. Da Rold and A. Pomarol, Phys. Rev. D 75, 055014 (2007) 
|arXiv:hep-ph/0612048| . 

Y. Hosotani and M. Mabe, Phys. Lett. B 615, 257 (2005) |arXiv:hep-ph/0503020| . 



A. Falkowski, Phys. Rev. D 75, 025017 (2007) |arXiv:hep-ph/0610336|. 

A. D. Medina, N. R. Shah and C. E. M. Wagner, Phys. Rev. D 76, 095010 (2007) 
|arXiv:0706. 12811 [hep-ph]] . 

M. Carena, A. D. Medina, B. Panes, N. R. Shah and C. E. M. Wagner, Phys. Rev. D 77, 
076003 (2008) |arXiv:0712.0 095 [hep-ph]]. 



K. Agashe, A. Belyaev, T. Krupovnickas, G. Perez and J. Virzi, arXiv:hep-ph/0612015 
B. Lillie, L. Randall and L. T. Wang, jarXiv:hep-ph/0701166| 
B. Lillie, J. Shu and T. M. P. Tait. larXiv:0706.3960l [hep-ph]. 



K. Agashe, G. Perez and A. Soni, Phys. Rev. D 75, 015002 (2007) |arXiv:hep-ph/0 606293 1 . 

S. J. Huber and Q. Shan, Phys. Lett. B 544, 295 (2002) | arXiv:hep-ph/0205327] . 

S. J. Huber and Q. Shafi, Phys. Lett. B 583, 293 (2004) | arXiv:hep-ph/0309252] . 

C. Csaki, A. Falkowski and A. Weiler, JHEP 0809, 008 (2008) [ arXiv:0804.l954l [hep-ph]]. 

S. Casagrande, F. Goertz, U. Haisch, M. Neubert and T. Pfoh, JHEP 0810, 094 (2008) 
|arXiv:0807.4"937l [hep-ph]]. 



[35] A. L. Fitzpatrick, L. Randall and G. Perez, Phys. Rev. Lett. 100, 171604 (2008). 



44 



[36] G. Perez and L. Randall, larXiv: 0805 .46521 [hep-ph]. 
[37] M. C. Chen and H. B. Yu. larXiv:0804.2503l [hep-ph]. 

[38] C. Csaki, C. Delaunay, C. Grojean and Y. Grossman, JHEP 0810, 055 (2008) [arXiv:0806.0"356l 
[hep-ph]]. 

[39] C. Csaki, A. Falkowski and A. Weiler. la7Xiv:0806.3757l [hep-ph]. 
[40] J. Santiago, JHEP 0812, 046 (2008) |arXiv: 0806 .12301 [hep-ph]]. 

[41] K. Agashe and G. Servant, Phys. Rev. Lett. 93, 231805 (2004) [ arXiv:hep-ph/0403143| . 

[42] K. Agashe, A. Falkowski, I. Low and G. Servant, JHEP 0804, 027 (2008) [arXiv:0712.2"455l 
[hep-ph]]. 

[43] G. Panico, E. Ponton, J. Santiago and M. Serone, Phys. Rev. D 77, 115012 (2008) 
|arXiv:0801.1645i [hep-ph]]. 

[44] E. W. Kolb and M. S. Turner, Front. Phys. 69, 1 (1990). 

[45] K. Griest and D. Seckel, Phys. Rev. D 43, 3191 (1991). 

[46] G. B. Gelmini, P. Gondolo and E. Roulet, Nucl. Phys. B 351, 623 (1991). 
M. Srednicki and R. Watkins, Phys. Lett. B 225, 140 (1989). 
M. Drees and M. Nojiri, Phys. Rev. D 48, 3483 (1993) |arXiv:hep-ph/9307208| . 
J. R. Ellis, A. Ferstl and K. A. Olive, Phys. Lett. B 481, 304 (2000) [arXiv:hep-ph/0001005] . 

[47] Z. Ahmed et al. [CDMS Collaboration], larXiv:0802.3530l [astro-ph], 

[48] J. Angle et al. [XENON Collaboration], Phys. Rev. Lett. 100, 021303 (2008) [arXiv:0706.0"039l 
[astro-ph]]. 

[49] |http://xenon.astro.columbia.edu/presentations/Aprile_IDM08.pdf| 
[50] T. Bruch [CDMS Collaboration], AIP Conf. Proc. 957, 193 (2007). 
[51] O. Adriani et al, arXiv:0810. 4995 [astro-ph]. 
[52] J. Chang et a/., Nature 456, 362 (2008). 

[53] See, for example, M. Cirelli, M. Kadastik, M. Raidal and A. Strumia. larXiv:0809.2409l [hep-ph] . 
[54] J. Hall and D. Hooper. lirXiv:0811.3362 l [astro-ph]. 



45 



